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Project  Description 


Quantum  transport  effects  including  electron  or  hole  tunneling  through  po¬ 
tential  barriers  and  buildup  in  quantum  wells  are  important  in  predicting 
the  performance  of  ultra-small  semiconductor  devices.  These  effects  can  be 
incorporated  into  the  hydrodynamic  description  of  charge  propagation  in  the 
semiconductor  device. 

Refs.  [1]  and  [2]  present  a  new  extension  of  the  classical  hydrodynamic 
model  to  include  quantum  transport  effects.  This  “smooth”  quantum  hydro- 
dynamic  (QHD)  model  is  derived  specifically  to  handle  in  a  mathematically 
rigorous  way  the  discontinuities  in  the  classical  potential  energy  which  occur 
at  heterojunction  barriers  in  quantum  semiconductor  devices.  The  model  is 
valid  to  all  orders  of  and  to  first  order  in  the  classical  potential  energy. 

The  QHD  equations  have  the  same  form  as  the  classical  hydrodynamic 
equations: 
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where  repeated  indices  are  summed  over  and  where  n  is  the  particle  density, 
u  is  the  velocity,  m  is  the  particle  mass,  Pij  is  the  stress  tensor,  V  is  the 
classical  potential  energy,  W  is  the  energy  density,  q  is  the  heat  flux,  and  Tq 
is  the  ambient  temperature.  Collision  effects  are  modeled  by  the  relaxation 
time  approximation,  with  momentum  and  energy  relaxation  times  Tp  and  r^,. 
Quantum  effects  enter  through  the  expression  for  the  stress  tensor  (and  for 
the  energy  density  derived  from  the  stress  tensor). 

Originally  the  quantum  correction  to  the  stress  tensor  in  the  QHD  equa¬ 
tions  was  given  to  0{Tt?)  and  involved  second  derivatives  of  the  classical 
potential. 

To  derive  the  new  effective  stress  tensor  and  energy  density,  we  construct 
an  effective  density  matrix  as  an  0{^V)  solution  to  the  Bloch  equation.  Then 
using  the  effective  density  matrix,  we  take  moments  of  the  quantum  Liouville 
equation  to  derive  the  QHD  equations  with  the  effective  stress  tensor  and 
energy  density  [Ij. 
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The  effective  density  matrix  has  the  form 
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where  V  is  given  in  center-of-mass  coordinates  R  =  |(x  +  y),  s  =  x  —  yby 
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Using  the  effective  density  matrix  in  the  moment  expansion  of  the  quan¬ 
tum  Liouville  equation,  we  obtain  the  QHD  conservation  laws  as  the  first 
three  moments  with 
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where  the  “quantum  potential”  is 
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The  quantum  correction  to  the  classical  stress  tensor  and  energy  density 
is  valid  to  all  orders  of  and  to  first  order  in  f35V,  and  involves  both  a 
smoothing  integration  of  the  classical  potential  over  space  and  an  averaging 
integration  over  temperature. 

We  define  the  ID  smooth  effective  potential  in  the  momentum  conserva¬ 
tion  equation  (2)  as  the  most  singular  part  of  U  —  Pn: 
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Figure  1:  Smooth  effective  potential  for  electrons  in  GaAs  for  50  A  wide 
unit  potential  double  barriers  and  50  A  wide  well  as  a  function  of  x  and 
logic  (T/300  K). 


The  double  integration  (over  both  space  and  temperature)  provides  sufficient 
smoothing  so  that  the  Pn  term  in  the  smooth  effective  potential  actually 
cancels  the  leading  singularity  in  the  classical  potential  at  a  barrier  (see 
Fig.  1),  leaving  a  residual  smooth  effective  potential  with  a  lower  potential 
height  in  the  barrier  region.  This  cancellation  and  smoothing  makes  the 
barriers  partially  transparent  to  the  particle  flow  and  provides  the  mechanism 
for  particle  tunneling  in  the  QHD  model.  Note  that  the  effective  barrier 
height  approaches  zero  as  T  — >  0.  This  effect  explains  in  fluid  dynamical 
terms  why  particle  tunneling  is  enhanced  at  low  temperatures.  As  T  oo, 
the  effective  potential  approaches  the  classical  double  barrier  potential  and 
quantum  effects  in  the  QHD  model  are  suppressed. 

Smooth  quantum  hydrodynamic  model  simulations  of  the  resonant  tun¬ 
neling  diode  were  presented  which  exhibit  enhanced  negative  differential  re¬ 
sistance  (NDR)  when  compared  to  simulations  using  the  original  0{J^)  QHD 
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model.  At  both  300  K  and  77  K,  the  smooth  QHD  simulations  predict  signif¬ 
icant  NDR  even  when  the  original  QHD  model  simulations  predict  no  NDR. 

The  ID  steady-state  smooth  QHD  equations  are  discretized  [3]  using  a 
conservative  upwind  method  adapted  from  computational  fluid  dynamics. 
The  discretized  equations  are  then  solved  by  a  damped  Newton  method. 

We  present  simulations  of  a  GaAs  resonant  tunneling  diode  with  A^Gai-xAs 
double  barriers  at  300  K  (77  K).  The  barrier  height  B  is  set  equal  to  0.1 
(0.05)  eV.  The  diode  consists  of  source  (at  the  left)  and  drain  (at  the 
right)  regions  with  the  doping  density  No  =  10^®  cm“®,  and  an  n  channel 
with  Nd  =  5  X  10^®  cm“®.  The  channel  is  200  (250)  A  long,  the  barriers  are 
25  (50)  A  wide,  and  the  quantum  well  between  the  barriers  is  50  A  wide. 
Note  that  the  device  has  50  A  spacers  between  the  barriers  and  the  contacts. 
We  have  chosen  parameters  to  highlight  differences  between  the  original  and 
smooth  QHD  models. 


Figure  2;  Smooth  effective  potentials  C/,  Ub,  and  Up  for  an  applied  voltage 
of  0.056  volts  for  0.1  eV  double  barriers  at  300  K.  x  is  in  microns. 
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Figure  3:  Current  density  in  kiloamps/cm^  vs.  voltage  for  the  resonant  tun¬ 
neling  diode  at  300  K.  The  solid  curve  is  the  smooth  QHD  computation  and 
the  dotted  line  is  the  0{ti^)  computation.  The  barrier  height  is  0.1  eV. 


Voltage 


Figure  4:  Current  density  in  kiloamps/cm^  vs.  voltage  for  the  resonant  tun¬ 
neling  diode  at  77  K.  The  solid  curve  is  the  smooth  QHD  computation  and 
the  dotted  line  is  the  0{h^)  computation.  The  barrier  height  is  0.05  eV. 
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Fig.  2  illustrates  the  smooth  effective  potentials  Ub  (for  the  barriers),  Up 
(the  Poisson  contribution),  and  U  (barrier  plus  Poisson  contributions)  for 
the  resonant  tunneling  diode  at  300  K  at  the  voltage  V  =  0.056  where  the 
I-V  curve  peaks  in  Fig.  3. 

Smooth  QHD  simulations  of  the  resonant  tunneling  diode  exhibit  en¬ 
hanced  negative  differential  resistance  when  compared  to  simulations  using 
the  original  0{h^)  QHD  model.  The  current-voltage  curve  for  the  resonant 
tunneling  diode  at  300  K  is  plotted  in  Fig.  3  and  at  77  K  is  plotted  in  Fig.  4. 
It  is  interesting  that  the  original  0{h^)  QHD  model  (see  Refs.  [4,  3]  and 
references  therein)  predict  very  different  I-V  curves — in  fact,  at  both  300  K 
and  77  K  the  original  0{h^)  QHD  model  fails  to  produce  negative  differential 
resistance  for  these  devices. 

In  his  lectures  on  Statistical  Mechanics  [5],  Feynman  derives  an  effective 
quantum  potential  by  a  Gaussian  smoothing  of  the  classical  potential.  After 
demonstrating  that  the  effective  free  energy  based  on  the  effective  poten¬ 
tial  is  accurate  for  smooth  classical  potentials  like  the  anharmonic  oscillator, 
he  goes  on  to  say  that  “it  fails  in  its  present  form  when  the  [classical]  po¬ 
tential  has  a  very  large  derivative  as  in  the  case  of  hard-sphere  interatomic 
potential” — or  for  potential  barriers  in  quantum  semiconductor  devices.  In 
this  investigation,  we  have  discussed  an  extension  of  Feynman’s  ideas  to  a 
smooth  effective  potential  for  the  quantum  hydrodynamic  model  that  is  valid 
for  the  technologically  important  case  of  potentials  with  discontinuities. 
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